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ABSTRACT 


An approximate theoretical analysis is made of a phase 
discriminator to be used as a synchronous detector of & very 
narrow-band FM wave in noise. Experimental work is reported 
which indicates the superiority of this method of detection 
over a conventional FM discriminator. 

A comparison is made of the effect upon limiter opera- 
ting point of two different methods of automatic gain control. 
A simple diode detector is shown to be superior to a "carrier- 
only" method when the carrier-to-noise ratio varies over a 
large range. 

The writer wishes to acknowledge the helpful assistance 
of the Study Group at Cubic Corporation and, in particular, 
that provided by Mr. Arthur FE. Noyes and Mr. Homer A. 
Lasitter. The continued assistance and encouragement of 
Professors C. F. Klamm, Jr. and P. E. Cooper of the U. S. 
Naval Postgraduate School is sincerely appreciated. The in- 
formal preparatory course which they conducted during the 
Fall Term constituted an indispensable background for this 


investigation. 
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l. Introduction. 

Cubic Corporation of San Diego, California, produces a 
Distance Measuring :quipment of the phase comparison type. 
When used for missile tracking the ground station transmits 
a carrier in the 200-400 Mc frequency range, frequency modu- 
lated by five lower frequencies. These frequencies are 
selected to permit unambiguous range determination out to 
one-half the wavelength of the lowest frequency. The mis- 
sile being tracked carries a transponder which receives the 
signal from the ground, extracts the modulating frequencies 
and remodulates a carrier of a different frequency which is 
transmitted to the ground. 

The accuracy of the entire system depends upon the 
phase shift experienced by the modulating wave as it passes 
through the transponder. A constant phase shift can be 
calibrated out by simply adjusting the range indicating 
servo at the ground station to a known range. Thus it is 
desirable to use an I.F. strip which exhibits a linear phase 
shift vs. frequency characteristic over a large range of in- 
put signal level. Van Voorhis has shown that synchronously- 
tuned stages are very tolerant to changes in tube capacity 
and their simplicity makes them desirable for this purpose 
(23)." Cubic's I.F. strips have six synchronously-tuned 
stages. The Q of plate tank circuits is deliberately made 
low (about 5.25) to decrease the steepness of the phase 


characteristic at midband and improve linearity at higher 


* Numbers in parentheses refer to bibliography. 








and lower frequencies; stated differently, the overall band- 
width of the I.F. strip is established by adjusting Q. 

To insure operation within the relatively linear por- 
tion of the phase vs. frequency characteristic, Cubic has 
considered it necessary to decrease the effective modulation 
index (6 و‎ to ae sufficiently small that only the 
first-order side-bands of the highest modulating frequency 
are significant. The method of negative feedback developed 
by Cheffee (4, 7, 8) is used to reduce B j.p, to approxi- 
mately..02. The feedback voltage is actually a portion of 
the re-transmitted signal which is offset from the incoming 
signal by the I.F. frequency of 30 Mc. 

The I.F. strip is followed by a limiter and Foster- 
Seeley discriminator. It has been Cubic's experience that 
the modulating signal power out of the discriminator must 
exceed noise power by about 10 db to permit accurate deter- 
mination of range. Since this deterioration of sigmal-to- 
noise ratio marks the limit of useful range, consideration 
was given to methods of improving performance. Several 
suggested changes were collected in an unpublished proposal 
submitted to the U. S. Air Force. Fig. (1-1) is a block 
diagram of the proposed transponder. It is the purpose of 
this paper to investigate certain of these suggested 


changes. 
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2. Characteristics of Random Noise. 

In recent years much has been written about the detec- 
tion of signals in noise. A rigorous treatment of any such 
problem requires the use of probability, since typical re- 
ceiver noise is a random processe 

This process cannot be written as sn explicit function 
of time, for even if its exact values over & long interval 
of time are known, future values cannot be predicted pre- 
cisely. It is possible, however, to determine the probabil- 
ity that the process will assume a certain value at any in- 
stant. This determination may be made from the observed 
values taken on by the process over an infinitely long in- 
terval of time, or by calculations when the characteristics 
of the primary source are known. Following the notation of 
Lawson and Uhlenbeck (16), this probability will be called 
the first probability distribution and is written W, (y,t), 
where y(t) is the random process and t is zs instant of 
time. The assumption is made that the E eek is station- 
ary, that is, any initial transients have decayed to negli- 
gibly small values and the process will display the same 
statistical properties over any (equal) long intervals of 
time. 

In a similar way an expression may be written for the 
probability that y(t) will lie in the range Y +A Y, at 
time ty and in the range Jot A y) at a later time i. 

This probability is called the second probability distribu- 


tion d is written W : 
end is written W(Y ,t.; y 


(7 or more simply 
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roy CT), where T = is t + Third and higher order 

probability distributions may be written, each describing 

the process in more detail. As Lawson and Uhlenbeck have 

pointed out on p. 35 of their book, this series of distri- 
butions completely defines the process. 

A convenient concept in this work is that of ensemble 
and time averages (16). If the random process y(t) is re- 
corded for a very long interval of time and then this re- 
cording cut into a large number (N) of intervals of length 
©, this ensemble of recordings may be considered the simul- 
taneous outputs of N identical receivers. Now the average 


value of a function of y at a time t may be written 


Oo‏ ہد 
Hp = dy f W (4%)‏ 
Oo |‏ — 
2709 ( 

where the bar indicates an ensemble average. This average 
me, in general, e function of t. If the average over time 
is performed the result is expressed as Ep. where the 
wavy line indicates a time average. If the averaging is 
carried out in the reverse order the result is expressed as 


E . The two methods yield the same answer since the in- 


finite time average of y(t) is obtained in either case. If 


the process is stationary, fos Fl) » Since 714) is 


Qe NUI Po 


independent of time and A (1) is the same for every member 


of the ensemble. 








Of special interest is the average value of y(t)y(t+7) 


Denoting this average by PIP), it may be written 
07 | 
AR) = y (©) (trt) = p p 
o 


is called the autocorrelation function of y(t), 


and for a stationary process becomes 
Oo 


br | | dy d^; 7۱ 42 We liza; v) 
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It is seen that p(o) = mean square value of y(t) 
D re 
If the average value of y is zero, &(0)= Q = variance. 


The importance of the autocorrelation function lies in 


the fact that He) and G(W), the power spectral density, 


are the Fourier r transforms of each other (249 
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In a large class of problems, the probability density 
functions for random noise voltages are relatively simple 
mathematical expressions. Lawson and Uhlenbeck (16) and 
Rice (19) are among those who have shown that typical re- 
ceiver noise sources, such as shot effect, thermal noise, 
induced grid noise, etc., are characterized by a normal or 


Gaussian probability density function of the form 
2 


W ( E | -2 Ha 
/ dv 9 JHT < o (2-4) 





Chapter 4 of Lawson and Uhlenbeck's book and Part II 
of Rice's paper contain extensive references to the litera- 
ture. Knudtzon (15) has published experimental verification 
for & variety of primary sources by measuring the autocorre- 
lation function using sampling techniques and then computing 
the power spectrum with an electronic differential analyzer. 

Two different representations of the random noise vol- 
tage will be found useful here. In the first form the vol- 
tage is called n(t), with a first probability distribution 


E" 
EL 
c. 


W (n) = E we (2-5) 
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In the second form the noise voltage is developed in a 
Fourier series occupying a narrow bend around the I.F. fre- 


quency £ Such that (TO) 


h(£)s X(2) coa 2NLÉ + Y(t) daca Tle , Q9 
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and ft, is measured from Lo. The a and b are independent 


random variables and x and y are orthogonal functions such 


that 








x (4) = uu uon 


Kid) y (trar zo for a & (2-7) 
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5. Effect of & Linear Device Upon Random Noise. 
It is seen from Eq. (2-3) above that if GW) is a 


constant for all values of c» , then 


ges Tas 


Dir) = ) 


a E = 
O ٣َ " ' ' کک٤‎ e m 


and the correlation between successive values of y is zero. 
This limiting case is called the "purely random process" 
(11) and the associated spectrum a "white spectrum". When 
such noise is modified by a linear device so that the output 
spectrum has the snape of the power frequency response of 
the device, the output voltage is still Gaussian (16) but 


the autocorrelation function is no longer discontinuous. 








4. Effect of Non-Linear Devices Upon Random Noise. 

Lawson and Uhlenbeck (16), Rice (19), Mideleton (17, 
18) and many others have published results of investigations 
of several types of non-linear devices under 5567 of 
signal and noise conditions. 

There are two generally accepted methods available for 
the study of the output of non-linear devices. In both of 
them the output power spectrum, at least in the range of 
interest, is determined. The first method, called by Rice 
the "direct method” and attributed by him to W. R. Bennett, 
Rn C. Williams, Je R. Ragazzini and others, is based upon 
the known theory of modulation in non-linear devices. In 
this method a typical modulation product is computed and 
the output power spectrum then obtained by summing all such 
products. 

The second method, attributed (19) to D. O. North, 

H. Van Vleck and D. Middleton, consists of first finding‏ .ل 
the autocorrelation function of the output, then taking the‏ 
Fourier transform to give the power spectrum of the output.‏ 
In this method use is frequently made of the characteristic‏ 
function to more simply evaluate the autocorrelation func-‏ 
tion. The characteristic function is widely used by Lawson‏ 
and Uhlenbeck throughout their book. Rice also uses the‏ 
characteristic function approach in many cases because of‏ 
its ready application to many devices expressible by the in-‏ 


tegral form 


10 





0 
dL. = P Le) E der , (4-1) 
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E 


where F(ju) is the Fourier transform of the I vs. V. charac- 
teristic of the device. 
In many problems one of the methods may prove simpler 


mathematically than the other. 
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5. Analysis of a Phase Discriminator. 


One of the changes suggested by Cubic is to replace the 
frequency discriminator with a phase discriminator which will 
be provided with a locally generated reference voltage in 
phase quadrature with the I.F. carrier wave. Cubic intends 
to supply this reference voltage by using an adaptation of the 
"active filter" designed by Jensen and McGeogh (14) to filter 
a portion of the I.F. signal which will then be shifted in 
phase by 90°. In this paper it will be assumed that the refer- 
Brice voltage is available. 

A sketch of the proposed phase discriminator is shown 


below in Fig. (5-1). 
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Fig. (5-1) 
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The voltages applied between the plates of V-l and V-2, 


respectively, and R.F. ground are 


C= Ep cot dat +E, are e re] + nt) 


Ez = Fr Cod tly £ ¬ Es hinf ev E ea] — Ge) 
or, since is small, 
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where n(t) has been expressed as in Eq. (2-6). 


Considering the two diodes to be ideal linear rectifiers, 


` the output voltage will be a high frequency wave with an en- 
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Terms after the third &re neglected since ihe condition 
will be imposed ٢ = 21 م‎ » It is recognized that values of 
IX] and ۳۹ equal to or greater than E. will occur, but 
within the range of our interest the probability of occur- 


ence will be very small. Hence 
=, = "۶ / + LX 


كيم - ما a‏ 


The binomial expansion is used again to expand the ex- 


pression within square brackets. 
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Now denoting by subscripts the quantities at E, and ور‎ 
BEES... x(t) x, y(t) = FL and letting X = و دلا اي تا‎ 
the autocorrelation function of the output voltage may be 


(5-4) 
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The noise in the I.F. amplifier işs taken to have a rectangu- 
lar power spectrum of width B 
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Since x and y are independent, their joint probability 


distribution function may be written as a product, e.g., 
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Calculation of the ensemble averares 4, Y uo ; 
and 3 9^ is most easily accomplished using the character- 


istic function method of Appendix (A). Thus 
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Before calculating d d. the assumption will be made 
that only one modulating frequency £ ) is present and that 
the effective modulation index is 3 Using the approxima- 
tion 
sin Y = Y= B un Zo) 7S 
the calculation is made simpler and a direct comparison may 


be made with the experimental results of Section (6). Thus, 
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The power spectrum of the output is found by taking the 
Using the results of 


Fourier transform of this expression 
Appendix (B), and writing the spectrum in positive frequen- 
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6. Comparison of Theoretical and Experimental Results for 
the Phase Discriminator. 


In the work at Cubic an accurate comparison was desired 
of the relative performance of the frequency discriminator in 
present use and the phase discriminator. Accordingly, labora- 
tory tests were conducted using an actual transponder I.P. 
strip and frequency discriminator. The latter was converted 
to a phase discriminator by disconnecting the end of C, marked 
S (see sketch below) from the primary of the transformer and 


feeding the reference voltage through C The reference volt- 


r 


age entered through & shielded cable, the lead of C, being 


1 
the only unshielded path in both circuit configurations. 
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Fig. (6-1) 


The complete laboratory set-up is shown in Fig. (6-2). The 
noise voltage n(t) was tube noise and thermal noise from the 
I.F. stages; carrier-to-noise (21) ratio (S,) was decreased 


by reducing the signal level. 
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For each setting of S» the output signal-to-noise ratio 


(8) was measured using the two discriminators in succession. 
All measurements were repeated on different days. Fig. (6-3) 


shows photorraphs of the oscilloscope presentations for vari- 
ous levels of 91" 
Eq. (5-8) may be used to determine $5.: In order to 


compare with the laboratory set-up, let 


+ = 9 6 

B = .27 

B= ¿Mc 
ZW - 6.6 Kc 


where 25 is the effective noise bandwidth of the post de- 


tection filter and is calculated from its measured 3 db 


. bandwidth using the table on p. 177 of Threshold Sigmals. 


This table also indicates that for the six single-tuned 


stages of the I.F. amplifier, the noise bandwidth closely 


approaches the 3 db bandwidth. 





Since S,= ھت‎ and f, 2.05 B, Eq. (5-8) may be used 
to give 


S =- SIGNAL Fe wen e 4, 
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Curves of > VS. S1 for R = 10 and R = 5 are computed 
from Eq. (6-1) and are shown in Fig. (6-4). If the addi- 


tional condition is imposed that 
7 A ee 
کک کی مه وس‎ 


then the asymptotic condition exists wherein 
S, = و‎ (6-2) 
where K depends upon P, . This curve is also shown in Fig. 
(6-4). 
Fig. (6-5) is a plot of experimental data showing the 
effect upon 84 of varying the ratio R. 

Fig. (6-6) is & comparison between theoretical and ex- 
perimental results when the conditions of Eq. (6-2) are met. 
Fig. (6-7) is a comparison of phase- and frequency- 

ı discriminator performance for a smaller modulation index, 
f, = .125. The theoretical result predicted by Eq. (6-2) is 


also shown. 
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7. Automatic Gain Control (ACC). 

The present transponder circuitry employs a limiter 
stage and AGC to the first four I.F. stages. The AGC detec- 
tor is an unbiased crystal diode. Both AGC and limiting are 
considered necessary because of the large range of input 
carrier level. The prosvect of increasing range even further 
by using the phase discriminator has aroused interest in an 


evaluation of different types of AGC. 


Of primary interest in this evaluation are the relative 
effects upon phase shift through the I.F. amplifier ana lim- 
iter as the carrier level varies over the range. The change 
in phase shift through these components is attributed to 
variations in the input capacity with bias voltage. A change 
in bias voltage alters the gain of the stage (which in turn 
changes the input capacity due to the Miller effect) and the 
direct plate current (which changes the grid-cathode imped- 
ance due to its effect on the size of the virtual cathode). 
Ine 1768691 1.۳ amplifier is considered sufficiently insen- 
sitive to these variations that it can be used with the 
present AGC or either of two proposed types of AGC. Accord- 
ingly, the manner in which limiter operating point varies 
with the strength and composition of the signal at its grid 
will be considered. 

The input circuit of the limiter stage is shown in Fig. 
(7-1). If the path from G to ground is regarded as an ideal 
linear rectifier in series with a fixed resistance Ry» an 


equivalent circuit may be drawn as in Fig. (7-2). 


E 








5702WA 


c‏ # مرس ل7 
(similar 10 OAKS)‏ رسد مب لس 
Ra UP‏ 
(3k)‏ 


Fig. (7-1)--Limiter input circuit 


10,4 
| M m 
b ۱ 


Fig. (7-2)--Equivalent circuit 


Since the grid path will conduct in only one direction 


end the capacitor will pass no direct current, any direct 


current flowing through the crystal must return to G through 


The d.c. grid bias is thus 


the ^K resistor. 
(7-1) 


The presence of Ex creates a biased linear rectifier. 


Let the input voltage be 
(7-2) 


S(L) + nz) = Ey tov wt + nt) 
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where the signal portion is considered unmodulated to make 
calculations more simple. Since the actual signal is fre- 
quency-modulated and a time average will be taken to obtain 
the d.c. component, the result obtained by using kq. (7-2) 
will be correct, 

The noise voltage is taken to be Gaussian with a first 


probability distribution 


Le 
۱ ^ AC* 
کد و‎ 13 
۷۷۰ qm ۶ 
Current flows through the diode when S(t)--n(t) - E, zO 


A= 
where E, is the d.c. bias voltage. Replacing ٢7۳۲7 


by the random variable Z and using the relation (9) 
۷۷)2( > ک ۳ ج۶‎ ) | ۱ (72a) 


where f(n) is the first probability distribution of n, the 
first probability distribution of Z is obtained, 


[226 - > 


Wia- ۱ E : 
(2) ۳ jan 3 $ (7-5) 


/ PA AAA MÀ کے‎ n Aun 
Now E E [309 * ^(£)-54] where the ensemble average 
c. x 


is taken over the region S(t)+n(t)-E,>0, or 


AAAA 


2 (7-6) 


_ | 
che 7 


over all positive Z. The ensemble average is taken first 
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The first integral on the right is integrated directly to 


give 
Ad ê DA ا‎ 
b. c. MECA 
2 6 = 
۳ ee 2G. = 
o 4/7 
- S(#) 


2 
Replacing x, by t? in the second integral on the right, 


and using the relation (12) 


7 7 
ce dr = ME 
2 


دن 


the result is 


MAA AAAA AAAA ANENA AN AN AN A 


B LEAT 
Ar INDIA حا ريس‎ sys [509-E] 


e- E 0 ےر‎ ELE p 2 B6-&] / ب سس‎ A 


ar 


5)4) 7 
Thu و8‎ VA ںار راے و وہ اہ‎ aA ج‎ 
ERA ہے‎ E 
T 07 
£ ور‎ E zm هوې‎ Zi y EL 2 €. E (7-7) 


Co 


From Eq. (7-1) and (7-6) a conditional equation can be 


written for Ep» 
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E, - (5t 269 MAAA AAAA AAAA 


eni (£)- 3‏ مهم هښ م٨۸۸‏ حصص مه هر٨‏ ۷ 


2 = 0ا 
gez e dt (7-8)‏ ا رر له 


It does not appear that Eq. (7-8) caen be expressed in closed 
form. Bennett (1), on p. 170 of nis paper, implied that nu- 
merical methods are necessary to solve this sel However, 
for a given value of carrier-to-noise ratio, the numerical in- 
tegration required to determine the time average is not too 
tedious. Values of I¿ o, for - 1000 >8 50 are of primary 
interest, since Middleton (17) and Rice (19) have written ex- 
pressions for, respectively, the no-signal case and very large 
carrier-to-noise ratios. In fact, with S(t) = 0, Eq. Gam) re- 
duces to Eq. 4-5, p. 785 of Middleton's paper. It will be 
shown later that for large S» results obtained from Eq. (7-8) 
are very nearly the same as those given by Rice on p. 125 of 
his paper (Vol. 24). 

An example will be used to illustrate the method of ne 
merical integration. Starting with Eq. (7-8) and choosing 


p 
E I و | ہا‎ the result is 


Eu سا‎ 
ANA AA AA 
Be BE j e R EM eet (7-9) 


IYWVAANAAA LAA INIA 
where use is made of the fact that dût =O. 
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In the actual circuit Ro is 3K<2. An arbitrary value of 


1K<L is assigned to Ry° This will affect the numerical re- 


sults but not the general conclusions to be reached. A value 





is now assigned to E and the two terms to be averaged 
are plotted on graph paper from wj t2- I to + JL in steps of 


2 2 
0.2 radians. The average value of each term is computed from 


the number of squares under the curve and the right side of 
the equation is thus reduced to a fraction of Ene The process 


is repeated for different values of EA 





until the right side 
is just E,. 

Values for Er were computed in this manner for carrier- 
to-noise ratios of zero, six, and ten decibels. In the tabu- 
lation below, E, is expressed as a fraction of P,, defined by 


E D 2 EZ 
zc OMS (7-10) 


and equal to the r.m.s. power which would be dissipated by the 


input voltage if it were impressed upon a one-ohm resistor: 


S- (db) Ep 
Less than -10 0552 Py 
0 .565 Py 
6 . 600 P. 
10 j^ 601 Py nen Integration) 
.607 2+ Rice's Approximation) 


The conditions for Rice's approximation are not fully met, 

since he required E,- Ej» O- and for 8j - 108b, حول‎ E, 230R 
An inspection of this tabulation reveals that the input 

voltage may vary from "pure carrier" to "pure noise" with an 


accompanying variation in limiter self-bias of less than 10 
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percent, provided that the mean square power input to the lim- 
iter is kept constant. 

Both of the proposed AGC systems are based upon circuitry 
which would present an extremely narrow bandwidth to noise. 
The system in Fig. (1-1) utilizes a md oscillator to 
translate downward in frequency a portion of the I.F. voltage. 
At this lower frequency (100 K,) the bandwidth may be reduced 
conveniently to 15 K,- The carrier-to-noise ratio after the 
bandwidth reduction (Sj) is then related to S, by 

2× / = E 
دور و هر‎ 5, 


or, 
A . 


S,/(yAs S,(db) - 2^2 db. 


Thus the AGC detector is only slightly affected by noise until 
the carrier is greatly reduced below noise. 

The other system is also a "carrier-only" type. Here the 
reference voltage رگ‎ conu, d would be shifted in phase by 909 
and multiplied by a portion of the I.F. voltage. The output 
circuit of the multiplier would have a very low cut-off fre- 
quency, hence a low equivalent noise bandwidth and the d.c. 
output would be proportional to the carrier level in the I.F. 
voltage wave. 

The simple diode which acts as the detector in the pres- 
ent AGC system has been investigated quite thoroughly by Rice 
(19). His Eq. (4.1-14) shows that for a quadratic detector 

Ty 2 <[o+ E] xa 


IC ۱ 
where X =constant. 


J9 








Using the notation of Eq. (7-10) this may be written 





Wa = < fz. (7-12) 
similarly for a linear rectifier 66 ( 1174 
TC 2 ZfirsyZ (5) +5 T(2) (7-13) 
dc. MAT ۱" 2 iet 
where I, and I. are Bessel functions of "و‎ argument. 
Again using Eq. (7-10) and the relation Spe E) Eq. (7-13) 


may be written 


La. GS 


= 


cao‏ اون esL)‏ ہے 


e 





Lm z| eÈ 
ee 
In Appendix C it is shown that the quantity within square 


brackets is very nearly unity over the entire range of Sy 


from 0 to CO . Hence, for the linear detector 
Li e CE RN 
If ideal AGC action is assumed, it is apparent that the 
crystal diode detector tends to maintain the total power out- 


put from the I.F. amplifier constant. The "carrier-only" AGC 


tends to keep Eg constant. From Eq. (7-10), by substituting 
2 
& 2س‎ the result is ': El 
= £s ( ۱ L) 
-~ سك ع | | کڈ‎ 
ee = 
Fa == 5). (7-15) 


For ود ان‎ decreasing S. requires that P, increase at a 


-2 
1 e 
the equipment will operate down in the region where 91 ap- 


l t 


rate proportional to S Cubic anticipates that eventually 


proaches -20db. It is apparent that in a system which has 
appreciable gain for large carrier-to-noise ratios, the re- 


quirement of Eq. (7-15) cannot be met with a practical 
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amplifier for Sy, << l. Jaffe and Rechtin (13) have pointed 
out this inapplicability of "carrier~-only" AGC to systems re- 
quiring a large dynamic range. 

There is another aspect of the AGC problem worthy of con- 
sideration. From Eq. (5-8) it is seen that the signal power 
out of the phase discriminator is proportional to e The 
phase discriminator is followed by a "compensating amplifier" 
of fixed gain which amplifies the five modulating frequencies 
before they are used to phase modulate a carrier. This car- 
rier is then multiplied up in frequency until it is 30 Mc 
below the incoming carrier frequency, & process DICE. 
multiplies the modulation index by the same factor (22). The 
modulation index out of the phase modulator at any given modu- 
lating frequency is proportional to the voltage at that fre- 
quency (3). Hence, the modulation index of the retransmitted 


wave depends upon E 5 . But 


8 = Pr - Pr (7-16) 


where Pr z modulation index at f4 of incoming signal. 
Br = modulation index at f, of retransmitted signal. 


Thus, if the circuit employs an AGC system which tends 
to maintain P, constant, a reduction in 91 will result in a 
smaller value of E,, which in turn tends to make 6 increase. 
Specifically, in the present system, > | , and for ڊ5‎ = 10, 
6 2 ۸2 ۰ If the phase modulator is assumed to exhibit a 
linear relation between modulating voltage and modulation in- 
dex, r= # E, 8 . Then from Eq. (7-16), 


/ 


PN 








Also, from Eq. (7-16), for the conditions given above 
Pr = Be-B,= #E (02) =. 98 1 
Thus, k= 49, where C is the value of Ea corresponding to 
C ^ 
و10 < و9‎ keeping in mind that P,= constant» 


Now from Eq. (7-15), solving for Eg 


“<= Mr 7) (7-18) 


= 2 P de m y 2ے‎ 
C سقف اسا‎ | = 1/2 py 4 


from which 


and, 











E 72 

AK, = ٢ س2207‎ 

Substituting this quantity for k, and using Eq. (7-18), Eq. 
(7-17) becomes 


8 eee (7-19) 
| hd y LL = $, S 
/ / 7476 BET د‎ S+] 
Evaluating Eq. (7-19) for 841- 1/100 or -20 decibels, 


= ج۰ا سلا‎ ES 
/ T S554 5 





١ 


From Eq. (6-1) it is noted that for large R, increasing 6, 

from .02 to .163 improves the signal-to-noise ratio S5 by 
2 

approximately (44e) — 66.5 or 87 


It remains to determine the effect of the I.F. amplifier 
phase shift characteristic upon the phase of the modulating 
wave at the output of the discriminator. Disregarding the 
noise voltage, the wave at the input to the first I.F. stage 


may be written (12) 
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CLA) = 4 رده‎ + @ ain ot] 
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Á |a dt + )S,(B)Aer[(oo+ no) Af 


> 
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19 
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DRAN sna ( 7-20) 


9 


کے 


The phase shift and amplitude response vs. frequency may 


be obtained from Eq. (5-8) p. 42 of Terman (22), 














(i+ SV ۱ 
(7-21) 
¡As | / 
7 “a. (248) 
As, ES 5 
where d K D 





E 
Then the wave at the output of the I.F. amplifier may be 


written 


e (4) = Ê | J (Ê air lst- e) 


gyenfeinajt= o + 4‏ د بط د + 


m 240 3 00.  ت ا زس‎ (7-22) 
nz 
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/ 
where ونا‎ Pp and ونا‎ P, are the amplitude responses and 
phase shifts at the frequencies fot rfi and fg7 nfi, respec- 
tively. 

Values of Jn (8) are tabulated below for a = .02 and 


.163 and n = A A11 higher orders are zero through at 
least the fifth decimal. 
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Values of by; ۸ ; br and 1, may be obtained from Eq. (7-21) 


and are tabulated below. 


بصت واف کھت A ape‏ سے ٦ nn‏ 9۰ ج۰ س 
سے — — nn ne‏ ہے دید سے e‏ ہے ھہ ہے oo mh deoa‏ ہو ہہس ا ہل ہہ ےد ے وے سے سے حص س 1 ارم mewn ere‏ سه در ده nn nenne Due‏ موس-ہےہ۔ جم‫ووے ھا ہے کے مس س e ta M‏ 
د سے ہہ ہے میسم —À——‏ —— سم وههمسه وې um nn -- ۰ u — — — —— —À‏ لوقصم جا ا یق ل ويم مر م عضا A SI‏ 








400° | 609 825 1 1 | سا مم 
پس مر سے ری O e‏ 
EUIS 852 — | 6‏ * ۲ 
uw ee er o mS Fei E ne pn men en c‏ ومون مه د مه م i- nd‏ 
ig | oO | 57. er 109.0 E‏ 
o E =‏ ۰ > — و 9 س 
n l 2 |‏ ۱ 





سے سس ات سان un‏ کو J— ——WLA——————————‏ ——————- ————————————— 


Taking b, = .02, it is seen from Eq. (7-22) that only the 
first order sidebands are significant in representing frequency 
modulation. <A phasor diagram of this situation is shown below 


‘in Fig. (7-3). 
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It can be seen from the diagram and Eq. (7-22) that the 


phase of the I.F. wave will pass through its zero position 


when the two first order sidebands are in such a position that 


their phasor sum has no component perpendicular to the don Wet 


phasor. 


ing 


‚343 ber کم۸)‎ +X) = 97 din (15° - x). 


from which 
din (Y 5 + 0 


J-to]. f Aer X‏ ہریرہ" 


a iS Cle x 


11 


linx = 


X 


0 
= 


This condition is expressed mathematically by writ- 


( 7-23) 


LI 7-7 > x) 
لا ور ےرس اٹ‎ 


A ری‎ 


sf ۸ 
و وہہ‎ 
ADS 
00193, 


/ 


Goo (7-24) 
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This is the amount by which the phase of the wave is shifted 
from the position it would have if the amplitude response were 
perfectly symmetrical about fg. Even if it were not negli- 
gible, it could be calibrated out. It is recorded here simply 
to illustrate the method of calculation and to serve as a com- 
parison with the phase shifts under different conditions. 
Again it is assumed that the I.F. amplifier is correctly 
tuned to fg and >, is increased to .163. The sidebands at 
n or 2f, must now be considered. A phasor diagram of this 


situation is shown below in Fig. (7-4). 





Fig. (7-4) 


It appears intuitively that the effective phase shift con- 
tributed by the unsymmetrical second order sidebands will be 
negligible because of their small amplitudes. To show that 


this is true, assume that the smaller one is removed 
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completely. The first order sidebands would have a nullify- 
ing component perpendicular to the lim uht reference for a 
value of x less than 10 of arc different from that of Eq. 
(7-24). Since the system accuracy as determined by the servo 
phasemeter is no better than 1.5 feet and 9 of arc corres- 
ponds to approximately 0.92 feet (at the ranging frequency of 
491 Kc), the shift in zero axis crossing due to the second 
order sidebands is obviously negligible. 

The above investigation assumed the I.F. amplifier to be 
correctly tuned to a center frequency fg. Now it will be 
assumed that the carrier frequency remains at Lo but that the 
center frequency of the I.F. amplifier is shifted upward to 
fot fq» This is an arbitrary choice but is based upon the 
fact that the 3db bandwidth is approximately 4f,, hence a 
shift in amplifier center frequency of fj narks about the 
maximum shift allowed for in the design "safety factor". The 
amplifier was of course designed for مم‎ z 402, so that only 
the first pair of sidebands had to be considered. 

Taking B = .02 and again disregarding the sidebands at 
fot 2f,, Eq. (7-22) may be written 


S.(0= B3.77 3,18) un [ut -el 


+ T (Brain fret -0 +60 
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Using the method of Eqs. (7-23) and (7-24) it is found that 


Xs 2060, (7-26) 

Now when £ = .163 it is seen that the amplitudes of the 
second order sidebands stand in almost exactly the same ratio 
as in the first case considered, while their phase angles 
differ by only 0.39. Thus, the effect upon zero crossings of 
the second order sidebands is again negligible. 

The amplitude mođulation which is present when ß = 1609 
is at twice the modulating frequency fj and is not passed by 


the tuned circuits of the compensation amplifier. 
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8. Conclusions. 

While the analysis of Section 5 is far from rigorous, it 
is considered a fair conclusion that use of the phase discrimi- 
nator in place of the frequency discriminator would provide a 
significant extension of system range. The experimental re- 
sults of Section 6 support this conclusion. In fact, the 
experimental results are likely to prove more typical in other 
discriminators of similar construction, since the assumptions 
and approximations involved in the work of Section 5 are 
really quite severe. The assumption that the diodes are per- 
fect switches, for example, is not realistic when the sigmal 
is very small. Linear operation of the diodes is similarly 
unlikely when the signal is very large. Complete symmetry of 
components is still another condition which is not attainable. 
However, it is believed that the mathematical model is a | 
reasonable approximation in the region where the r.m.s. refer- 
ence voltage is maintained very large compared with the total 
r.m.s. voltage in the I.F. amplifier [conditions of Eq. (6-2) 
or carrier-to-noise ratio (84) is large. This condition is 
quite easily met when total r.m.s. voltage in the I.F. ampli- 
flier is maintained nearly constant. The reference voltage 
emplitude is simply kept constant at a value which is large 
with respect to the carrier amplitude where $4 =+10 db. This 
is a strong point in favor of using the "constant total power" 
AGC method of Section 7. 

The desirability of keeping circuit changes to a minimum 


suggests another advantage of retaining the present AGC 


45 





system rather than using a "carrier-only" type. Even the ad- 
dition of an amplifier for the AGC voltage would be relatively 


simple compared with the extra circuitry required, for ex- 


ample, in Fig. (1-1) ۰ 


46 








9. Recommendations for Further Work. 


It appears that the phase discriminator will find many 
applications in this send other areas of research where phase 
measurements and/or correlation detection are employed. Radio 
navigation, radar, color television and many other electronic 
devices could be cited as examples. Consequently, a really 
thorough investigation of the phase discriminator might be a 
highly useful project. 

It was mentioned in Section 4 that many non-linear de- 
vices have been investigated theoretically to determine their 
effects on the spectrum of sigmal plus noise. Regardless of 
the method used, a fairly common occurrence is the appearance 
of an integral or infinite series which cannot be expressed 
in closed form. 

Blachman (2), Middleton (17), Rice (19) and Lawson and 
Uhlenbeck (16) are a few who have reported such results. 

Their procedure is then to assume "very large carrier-to-noise 
ratios" and "very small carrier-to-noise ratios" which reduce 
the equations to closed forms. These results are useful for 
some practical applications such as the detection of extremely 
weak signals in noise or circuit response to signal alone. 
However, equally valid results are often obtainable by the 
methods of modulation theory for these limiting cases. The 
region in which calculations are most difficult (carrier-to- 
noise ratios roughly -20 db to +20 db) would therefore be the 
region of most practical interest. 

It would then seem that a logical approach to a thorough 


investigation of the phase discriminator problem would involve 
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a determination of an integral or series form most amenable 

to solution on a digital computer. Two series might be re- 
quired to cover the entire range (17) but as long as they 
could be derived from the initial model without approximations 
the numerical results could be obtained to any degree of ac- 
curacy desired. 

Eq. (4-1) requires the knowledge of F(ju) for the device 
being considered. Rice (19) has tabulated F(ju) for a number 
of non-linear devices in Appendix 4A of his paper, but no one 
has written this function for the phase discriminator. This 
function could probably be calculated with the aid of the 
papers by Rice (19), Rice and Bennett (20) and Middleton (17). 
Having obtained this F(ju), the simplest form for numerical 


analysis would have to be reasoned out. 
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APPENDIX A 
A.1 AVERAGE VALUES 


The characteristic function for & random variable y, or 


for the corresponding distribution, is defined by (10) 


ve; ; TA 
(= e = [dee Fl ee 


Eos 


This result may be generalized for a distribution of n 
variables (10) but for the present purpose it is sufficient 


to write 
Zz 
nen Ps EX Y) dy dy. گت‎ 


Since ty and t, are constants as far as the integration 
۱ concerned, it ig possible to differentiate both sides with 


respect to با‎ and Vr Differentiating m times with t, and n 


times with E the result is 


E af ( W, (a 2ل‎ EXP 
5 E E | ٠ 
/J (£, "Zn tay a) d» d, : 


If now uu and t, are set a to zero, the result is 


EU رل‎ (4, 42) men by >) olu of 
E | fretis ER 
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MNFN —— 
The right side of this equation is just (4) u” Y, 
/ e 


Hence, 


epp 


0 f i (1) 


sd H(t 4) 
d 2£ ۵ f, d 





(A.1-3) 


٨د‎ ٥ : 


A.2 CHARACTERISTIC FUNCTION FOR Wo 


The characteristic function for W, will now be calcula- 


ted. 

TETEA = | | lara i eg 2 
tee zc + 4 -2pg 2 3 

= بې‎ Je y plian- zat] 


E 2 
JA expl 14, ۲ ze پا‎ -2 2 


The integration over Yi is performed first. Completing the 


square in the exponent the result is 


D E / 2 2 2 22 
d Expl J 2 A 7 2 5 uix Stes 
LE ۳ ۱ / 8 3 

ل alt‏ ےو ۔ رت 7 EC‏ 
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Now making the substitution 


و6 - رل ع< 22 ١‏ سے 


this integral reduces to i 


2 
reg exp/ IRR ge LIE 
2 7۰۰] مج‎ tA eg 5 م2‎ 
OO 


Making use of pair No. 708.0 Campbell and Foster (6) this be- 








E 2) e EXP PB 4- رج‎ 58 re 
| 1 un ۱ 2 
end 
y te 
VG, t,)= اس‎ M expl. 1 EIL 
1 6 2 - 
۳ د‎ P> +4 Apr - d 
Log) ته‎ 1 
0 سور‎ — #2, 
Ma (ee 0 
- 00 
Setting 
322-22 0-2 


end again using pair No. 708.0, this becomes 


e. z a 2 2 
pE t= EZ ELA TE (T+A) 








which simplifies to 
- £[£^-2e2. 2, 27] 
y CE A) = € (A. 2000) 
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APPENDIX B 


Transformation of Eq. (5-7) requires the determination 


of the Fourier transforms of Cox. a), ^ 2 Cox CU, 9 
۱ 2 


p ceo e 3 l p* gnd pP” 
As a starting point, it is known that (16) 
CO N 
AZUFL 
لم‎ E = 0)۶ (B-1) 


- 00 
where Ó(f.) is the Dirac delta function and has the following 


properties: 


Yo سرت‎ 


ہج ت77 CO,‏ 


"n ff = |. 
SID = IA 


۲ Jr ۳ 


To transform Cd AT, it is first written 
-۷ 5 
بھی‎ yv m L| et رہ‎ ۱ 
Ce 
50 eo ۰ زرم‎ 
E D - u^ 
F [cover] 2 | له عص‎ yr 
~ O0 


ES ceo 
2 TE e S -IRTEE 
- je و‎ E تو‎ 04070 


OO CO 
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"E 6(-5* £ Sed) 


or 
F [cea uwr] = á4 - f). (positive frequencies) (B-2) 


The transform of eq) is obviously 


Mea | (Ma, ٣٣ 
PECES (B-3) 
oe elsewhere d 


The transform of PX) is most easily taken by using 


the convolution theorem (25). Thus, 


J 808 


Gl) G(f-2) d4 


— 
qud 


(B-4) 


اد 


[ EO <O. 


DA 
— | Pu 


Ol 
m 


Denoting the transform of PT) by H(F) and using 


the convolution theorem again, the transform of 249 is 


found to be 
co 
-J2TAT / | 
Je € حر‎ per HO Gu (f -4) da 
en Pp 
2 ٦ 
Ep. -— 7 9 >> 2 8 
عم‎ 8 2 m 
O eliewbhsre, (B-5) 
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It was seen above (lq. B-2) that transforming Cod wv 
shifted the spectrum of the unit impulse by tf . In an 


entirely similar manner it is found that 


= E y o«f«B -f 
Aut 
"نه مه م‎ ۱ d = = ۱ E -f<f< Ë +f (8-6) 
290 Oa elsewhere. 
p 2f, OC sE <£ 
2 
co 
“mtr 2-3 f£,« £« &f 
0 B) 
لے‎ de dr = (B-7) 
pj £^, B-f« £4 Bf 
p! | 
260 
Ó ) elsewhere. 
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APPENDIX C 


-X 
Values of 1,(x) and 1,(x) for 05X5 Sand of € IS 
-X 
and € کو‎ for 6£ x €20 may be obtained from British 
Association for the Advancement of Science Mathematical 


Tables, Volume VI. (5). 


The expression 
S, 


7/3) 5 — —— (1-5) T,(5) e رک‎ 


was evaluated in 3 db increments of S- from -15 db to +15 


db. Results are tabulated below: 





5, (db) Sy T 
-15 02125 1.000 
E .0625 1.000 

x .125 1.000 
-6 .25 1.000 
-3 .5 1.008 
0 1 1.024 
5 2 1.045 
6 4 1.073 
9 8 1.099 
12 16 - 1.110 
15 32 1120 


The asymptotic values of Ij (x) and I (x) are found on 


۳ 2/1 of the same tables. Substituting these expressions 


into the right side of the equation for a? and letting 


DID 





>= x , the result is 


a 
T(S)) = نك‎ [Ir ا وک‎ 
) ۷ ATX Bx t 2)7 


(a ک‎ NS 6 my os 2 سو‎ n 7 


Now letting لز‎ —> wand evaluating the limit, 


: SS : ER مز و‎ 
En م8‎ ) Lin | Gun 
رک‎ oo 
E 4 = CIENT 
x=>0 ZT 2 


We TRUST 
X 
E | / 


Zu 


meat û‏ مسیون موس 


له 
1 








|| 
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For ©2,=0O 


TO = as : 
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